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Abstract 

| Some problems of the perturbative quantum chromo dynamics (pQCD) 

jet calculus are discussed. The first one is related to the terminology of 
the order of calculation. Due to cancelation of LO and NLO terms in 
the ratio of mean multiplicities in gluon and quark jets r the nowadays 
obtained results about it should be called as 4NLO approximation. 
The second problem reveals itself in calculations where corrections to 
some values (in particular, to r') are larger at present energies than 
lower order terms. Some characteristics which do not suffer from this 
■ deficiency are proposed. Next problem lies in interpretation of the 

OO \ negative values of cumulant moments which are considered as an in- 

dication to the replacement of attraction by repulsion in sets with 
definite particle contents. Finally, the problem of the generalization 
, of QCD equations for generating functions is briefly discussed. 

o 

r-| ■ The numerous achievements of pQCD in prediction and description of 

Ph. properties of quark and gluon jets are well known and described in many 

P-i! review papers (see, e.g., |], 0, |j, |5|). Here, I would like to discuss some 

problems, related to these calculations and, often, left behind the scene. 

First, let me remind some simplest definitions |], |I[ concerning jet mul- 
tiplicities in QCD. The generating function G is defined by the formula 



X 



G(y,«) = EW. C 1 ) 



n=0 



where P n (y) is the multiplicity distribution at the scale y = \n(pQ/Q ) = 
ln(2Q/Q ), p is the initial momentum, B is the angle of the divergence of 
the jet (jet opening angle), assumed here to be fixed, Q is the jet virtuality, 
Qo = const, u is an auxiliary variable. 

The moments of the distribution are defined as 



p = E n Pnn(n-l)...(n-q + l) _ 1 d*G{y,u] 
q (EnP n n)i (n)*' du" 



u=l ) 



K 1 dnnG(y,u) 
9 (n)' dui 



u=l- 



(2) 



(3) 



Here, F q are the factorial moments, and K q are the cumulant moments, 
responsible for total and genuine (irreducible to lower ranks) correlations, 
correspondingly. These moments are not independent. They are connected 
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by definite relations which can easily be derived from moments definitions in 
terms of the generating function: 

F q = £ C™_,K q _ m F m . (4) 

m=0 

The QCD equations for the generating functions are: 

G' G = f 1 dxKg(x)>y*[G G (y + kix)G G (y + ln(l - x)) - G G (y)] 
Jo 

+ n f f 1 dxK^x) 1 2 [G F (y + \nx)G F (y + \n(l~x))-G G (y)], (5) 

G' F = f 1 dxK<£( x )^lG G (y + lnx)G F (y + ln(l - x)) - G F (y)\, (6) 
Jo 

where G\y) = dG/dy, rif is the number of active flavours, 

7 = ——i ( 7 ) 

71 

the running coupling constant in the two-loop approximation is 
where 

_ llN c -2n f 17N 2 -n f (5N c + 3C F ) 
Po — ^ , Pi — ^ , (?) 

the labels G and F correspond to gluons and quarks, and the kernels of the 
equations are 

K£(x) = --(l-x)[2-x(l-x)], (10) 

Ob 

K F G {x) = ^[x 2 + (1 - x)% (11) 



--1+- 

x 2 



(12) 



iV c =3 is the number of colours, and C F = (N 2 - 1)/2N C = 4/3 in QCD. 

Herefrom, one can get equations for any moment of the multiplicity dis- 
tribution both in quark and gluon jets. One should just equate the terms 
with the same powers of u in both sides of the equations. In particular, the 
equations for average multiplicities read 

(na(y))' = J dx-f 2 [KS(x)((n G (y + lux)) + (n G (y + ln(l - x)) - (n G (y))) 

+n f K£(x)((n F (y + In a;)) + (n F (y + ln(l - x)) - (n G (j/)))],(13) 
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(My))' = J dx^K^{x){(n G {y + lnx)) + (n F {y + ln(l-x))-(n F (y))). (14) 

Their solutions can be looked for as 

ry 

(n GjF ) oc exp(y 7 G Ay f ) d y f )- ( 15 ) 

Using the perturbative expansion 

1g = 7 = 7o(l - ai7o - fl27o ~ «s7o) + O(7o)> ( 16 ) 
one gets the solution in the form || |7|, |j 

(n GtF ) = A GtF y~ a ^ exp(2 Cv ^ + S GiF (y)), (17) 

where c= (4iV c //3b) 1/2 , 



-3/2\ 



* Q (v) = ^[2a 2 c 2 + %(\n2y + 2)] + -[a 3 c 2 - ^(ln2y + 1)]) + 0(y 

y/y Po y Po 

(18) 

Usually, in place of j F the ratio of average multiplicities in gluon and quark 
jets 

r = ^ = ^exp(fo(y)-My)) (19) 

is introduced, and its perturbative expansion 

r = r„(l - r l7o - r 27o 2 - r 3 7 3 ) + O( 7o 4 ) (20) 

is used. The analytic expressions and numerical values of the parameters 
di, rj for alH < 3 have been calculated from the perturbative solutions of the 
above equations (the review is given in 0). Within these approximations 
the experimental data about mean multiplicity in e + e~-annihilation are well 
described as seen in Fig. 1 where the notation K = 2A F = 2A G /r Q is 
used.. However the data about the ratio r can be described with much lower 
accuracy about 15% in such an analytical approach (see Fig. 2) even though 
each subsequent perturbative approximation improves the agreement. 

However, one should mention here the quantitative fit provided by the 
computer solution of the equation || |K|. This poses the question about the 



accuracy of perturbative approximations for this particular characteristics 
and indicates that the higher order corrections are still comparatively large 
for this ratio up to the highest presently available energies. Let us also note 
that the exact solutions of these equations for fixed coupling constant were 



given in |TT], 

The relation between the anomalous dimensions 7 of gluon and quark jets 

is 

If = 7 - -, (21) 
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where 

r > = dr/dy = Sr r l7o 3 [l + ^ 7 o + (7^ + #i)7o + O( 7o 3 )] (22) 

with r = § = 9/4; 5 = /3 /8iV c ; = A/4iV c /5 . 
Thus 

7F = 7o [l-a l7o -(a 2 +Br 1 )7 2 -(a3+2Br 2 +Br2) 7o 3 -(a4+B(3r3+3r 2 r 1 +B 1 r 1 +r?))7^]. 

(23) 

In these expressions we meet with two problems. 

• Terminology. 

The two leading terms in the energy behaviour of quark and gluon jet 
multiplicities are absolutely the same as seen from Eq. ([T7|) and cancel 
in their ratio r fll9|). Therefore this ratio is given by tq = 9/4 both 
in the leading (LO) and next-to-leading (NLO) approximations. Thus, 
the common notation DLA, which is used in Fig. 2 near the value 
r=9/4, should be considered as LO+NLO-prediction of QCD for the 
ratio r. Therefore, the term r l7o in (p0[ ) describes 2NLO corrections to 
the anomalous dimension. However, in the literature, it is often called 
as a MLLA (NLO) term what is wrong. Nevertheless, namely such 
notation is commonly used in Figures. Here, in Fig. 2 we have used 
the notation with the letter r added at the end. It implies that, e.g., 
3NLOr means that the term with 7 g in the perturbative expansion of r 
has been taken into account but it corresponds to 4NLO-contribution 
to the anomalous dimension. 

A misuse of the terminology for the anomalous dimensions 7's and for 
the ratio r is clearly displayed in the explicit expression for 7^ (|23|). 
Its last 4NLO term contains 04 which has not yet been calculated. 
Together with it, the contribution from r is present with all terms cal- 
culated already and containing for i < 3 only. Thus, let us stress 
again that in this sense one should say that "r 3 "-term in r corresponds 
to 4NLO contribution to the anomalous dimension of the quark jet even 
though it is proportional to 7o in the perturbative expansion of r. 

• Calculations. 

The cancellation of two leading terms in the ratio r reveals itself also in 
the proportionality of the scale (energy) derivative r' to 7$. Therefore 
it can be calculated up to the terms 0{pfl). The leading term is very 
small (about 0.02 at the Z°-resonance). Asymptotically, all corrections 
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vanish. However, at present energies of Z°, they are still quite impor- 



tant. The second term in the brackets in (22) is larger than 1 since 
2r 2 /ri « 4.9 and 7 ~ 0.45 — 0.5. Even the third term is approximately 
about 0.4. The problem of convergence of the series at Z°-energies and 
below becomes crucial. The derivative of the ratio r (its energy slope) 
is very sensitive to high order perturbative corrections. 

Therefore, it is desirable to use at present energies such characteris- 
tics which are less sensitive to these corrections. In particular, these 
corrections partially cancel in the ratio of derivatives (slopes) 

r (D = m:. ( 24) 

(n F )' 1 ; 

The same is true for the ratio of curvatures (or second derivatives) 

.( 2 ) _ m: 



(25) 



(n F )"' 

The QCD predictions for them 

r < r {1) < r (2) < 2.25 (26) 
were recently confirmed in experiment (see Figs. 3, 4 from [IB]). 



Interpretation. 

Another question I'd like to raise concerns physical interpretation of 
oscillations of cumulant moments in QCD which is not yet completely 
clarified. Usually exploited phenomenological distributions of the prob- 
ability theory do not possess any oscillations. E.g., all cumulant mo- 
ments of the Poisson distribution are identically zero. One interprets 
this as the absence of genuine correlations irreducible to the lower-rank 
correlations. For the negative binomial distribution one easily gets 

K 2 

Hq = til = > 0. (27) 

" F q q(q + l) V } 

Since F q are always positive according to their definition, this inequality 
implies the positive values of K q . 

In the leading order approximation, the gluodynamics equation for the 
generating function 

[logG(y)r = 7o 2 (G(y)-l) (28) 



5 



transforms in the relation 



q 2 K q = F q or H q = ^. (29) 
However already in the next-to-leading order if g -moments become neg- 



ative with a minimum at the rank q min fa + 0.5 ~ 5 [|TJ[]. This 
minimum is rather stable. It slowly moves to higher ranks with en- 
ergy increase and disappears in asymptotics as is required according to 
the formula (p9|). At higher orders of the perturbative expansion, the 
oscillations of higher rank cumulant moments show up They are 
confirmed in experiment |16|, [17|] (see Fig. 5). 

Let me mention here that the plots of D q = q 2 H q instead of H q would 
be even more instructive to reveal the oscillations. In this case they can 
be easily compared to the LO prediction according to which = 1. 
Also the comparison to results of the negative binomial distribution 
would become simplified. The plot of NBD results shows monotonic 
increase of D^ BD from 1 at q = 1 to 2 at q — > oo which is significantly 
different from QCD oscillations. 

Both the role of conservation laws and the changing character of the 
genuine correlations can be blamed as originating these oscillations. If 
the latter factor is important it would imply that attraction (clustering) 
is replaced by repulsion (and vice versa) in particle systems with differ- 
ent number of particles. It would be interesting to find other examples 
of such a behaviour in hadronic systems. 

Generalizat ion . 

Finally, there exists the problem of possible generalization of the equa- 
tions for the generating functions. From one side, we understand that 
even if treated as kinetic equations these equations are limited by our 
ignorance of non-perturbative effects, simplified treatment of conserva- 
tion laws etc. Some phenomenological attempts to avoid these limita- 
tions were attempted from the very beginning [19|, |20|. In Jl8| it 



was proposed to treat hadronization of partons at the final stage of jet 
evolution in analogy with the ionization in electromagnetic cascades 
where it leads to their saturation and to the finite length of the shower. 
Three different stages of the cascade were considered in the modified 



kinetic equations proposed in [19], pOfl . No quantitative results were, 
however, obtained. 

The most successful modification of above equations was recently pro- 
posed in the framework of the dipole approach to QCD with more 



accurate kinematic bounds. It has been shown that the ratio r can be 
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obtained in good agreement with experimental data. Nevertheless, fur- 
ther study [22[ of higher rank moments of the multiplicity distribution 
predicted by the modified equations has shown their extremely high 
sensitivity to higher orders of the perturbative expansion. As shown 
in Fig. 6, the moments diverge at high orders and the only trace of 
oscillations can be noticed in the changing signs of the moments of 
the subsequent ranks. The results become inconclusive. Thus no suc- 
cessful generalization is at work nowadays. Rather, the general trend 
shifted to the direct calculation of non-perturbative effects in some jet 



characteristics (see, e.g., p3 



At the same time, the success of numerical solutions of the existing 
equations [0, 10 1 raises the question if the generalization will give any 



other noticable contribution and our failure to describe more precisely 
the ratio r could be just due some defects of the purely perturbative 
expansion at available energies. More rigorous treatment of the nu- 
merical solutions of the equations should be done. Moreover, it was 
claimed recently 03 that the renormalization group improvement of 
the perturbative results gives rise to good description of experimental 
data. 



In conclusion, I'd say that, even though some principal questions concern- 
ing the calculation of some properties of quark-gluon jets and the validity of 
QCD equations for the generating functions at higher orders are not yet 
resolved, the practical accuracy of the pQCD calculations is high enough, 
especially, in view of the rather large expansion parameter. 



This work is supported by the RFBR grant 00-02-16101. 



Figure captions. 



Fig. 1. The energy dependence of average multiplicity of charged parti- 
cles in e + e~-annihilation. The results of different fits according to formulas 
of perturbative QCD and of the Monte Carlo models are shown ( the solid 
and dotted lines are the fits of formula (|18D with one and two adjusted pa- 
rameters, the dashed line is given by the HERWIG Monte Carlo model; the 
vertically shaded area indicates the gluon jet data multiplied by the theoret- 
ical value of the ratio r ( p0|) ) . 



Fig. 2. The experimentally measured ratio r of multiplicities in gluon 
and quark jets as a function of energy in comparison with the predictions of 
analytical QCD and of the Monte Carlo model HERWIG (different QCD ap- 
proximations, described in this paper, as well as r(e) with integration limits 
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e y and 1-e y in Eqns (Q), (||) are indicated at the corresponding lines). 

Fig. 3. The ratio of the slopes of the energy dependences of mean mul- 
tiplicities in gluon and quark jets according to experimental data and some 
theoretical calculations. 

Fig. 4. The ratio of the curvatures of the energy dependences of mean 
multiplicities in gluon and quark jets according to experimental data and 
some theoretical calculations. 

Fig. 5. The measured ratio H q of the cumulant and factorial moments 
oscillates as a function of the rank q according to experimental data on 
multiplicity distributions of charged particles in e + e~-annihilation at the Z° 
energy (the inset in the upper right corner shows the data for the moments 
of the ranks 2, 3 and 4). 



Fig. 6. The if g -moments in the modified dipole approach |2T, |22| dras- 
tically diverge at higher orders for large ranks q with changing the sign at 
subsequent ranks. 
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